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Abstract 

The universal formulation of spin exchange models related to Calogero-Moser mod- 
els implies the existence of integrable hierarchies, which have not been explored. We 
show the general structures and features of the spin exchange model hierarchies by 
taking as examples the well-known Heisenberg spin chain with the nearest neighbour 
interactions. The energy spectra of the second member of the hierarchy belonging to 
the models based on the A r root systems (r = 3, 4, 5) are explicitly and exactly evalu- 
ated. They show many many interesting features and in particular, much higher degree 
of degeneracy than the original Heisenberg model, as expected from the integrability. 



In a previous paper [[]]] we presented universal Lax pairs for the spin Calogero-Moser 
models [Q, |3], |4j and spin exchange models || |6|, 0, |] which give rise to enough number 



of conserved quantities. As in the continuous Calogero-Moser models [pi fLO, O, these 



conserved quantities are shared by all members of the integrable models belonging to the 
same hierarchy. However, in the context of spin exchange models which we will mainly 
discuss in this paper, the notion of integrable hierarchies is rather unfamiliar and virtually 
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unexplored in contradistinction with the soliton hierarchies or Calogero-Moser hierarchies. 
The purpose of the present paper is to call attention to this interesting and potentially 
useful subject in mathematical/statistical physics by examining the second member of the 
integrable hierarchy of the simple and well-known Heisenberg spin chain with the nearest 
neighbour interactions. The energy eigenvalues of these models are evaluated explicitly and 
exactly for relatively small numbers of the spins and remarkably high degree of degeneracies 
of the eigenvalues is observed. We believe that this fact is closely related to the existence of 
a large number of conserved quantities shared by the hierarchy. 

The hierarchy of the spin exchange models consists of models associated with a fixed 
root system A and the totality of irreducible representations of the finite reflection (or 
Coxeter or Weyl) group In the case of the Heisenberg spin chain the root system is 
A r , which corresponds to the Lie algebra of su(r + 1). In this case the Weyl group of A r is 
the well-known symmetric group S r+ i. First let us rephrase various concepts and dynamical 
quantities of the well-known Heisenberg spin chain in terms of the languages of the weights 
and roots and the associated reflections, which are useful for the universal description of 



the Calogero-Moser models [12], [H| |T|, [16| and the spin exchange models |l|]. If the A r 
root system is embedded in R r+1 , it has the following simple representation in terms of an 
orthonormal basis {e,,} of R r+1 : 

A r = {e j -e k \j,k = l,...,r + l}. (1) 

The original Heisenberg model is based on the vector (V) representation of A r corresponding 
to the simplest Young diagram fj. It consists of the following weights: 

□ : V = {e J+7 | j = l,...,r + l}, (2) 

in which 7 is orthogonal to all the A r roots and it will be omitted hereafter. At each "site" 
ej (or simply the site j) a spin space is attached. Throughout this paper we consider the 
simplest case of the 1/2 spins of su(2), that is ups (|) and downs (|). In other words, 
the state space of the original Heisenberg chain is a 2 r+1 dimensional vector space. The 
Hamiltonian now reads 

H = JJ2V P , JeR, (3) 

penw 

in which J is the coupling constant and n^ 1 ) is the set of affine simple roots of A r : 

n (1) = {ei - e 2 , e 2 - e 3 , . . . , e r - e r+1 , e r+1 - ej, (4) 



and V p is the spin exchange operator corresponding to the root p. In the present case P e ._ e . +1 
exchanges spins at the sites j and j + 1 (to be more precise the sites ej and e J+1 ) and keeps 
the other spins intact. Let us denote the Pauli spin matrix acting on the site simply by 
<7j, then V p has a well-known realisation: 

p = ej - e k =>■ V p = (1 + <Xj ■ o ! fc)/2. (5) 

Thus we find that the above Hamiltonian describes the nearest neighbour Heisenberg 
spin chain and the last root in (||) provides a periodic lattice. 

It is worth remarking that the above Hamiltonian (H) can be obtained as a special limit of 
a more general spin exchange model with elliptic exchange functions |§, and the conserved 
quantities can be constructed |§ as shown by one of the authors. The above Hamiltonian (|3|) 
and its conserved quantities are shared by all the members of the same integrable hierarchy. 
Each member of the hierarchy, to be called a model TZ for short, is specified by an irreducible 
representation TZ of the Weyl group of A r . The elements of TZ, to be called "sites", too, are 
vectors in R r+1 which form a single Weyl orbit: 

P e TZ =>• s p {p) eTZ, Vp g A, (6) 

in which s p reflects a vector £ G R r+1 in a hyperplane perpendicular to the root p: 

^s p (0 = Z-(p v -OPi P y ^lplp\ peA. (7) 

As above, an su(2) spin, an up (f) or a down (!) is attached to each "site" p. Thus the 
state space of the model TZ is a 2 D dimensional vector space in which D is the total number 
of sites in TZ. The action of the spin exchange operator V p is now 

Vp : spin at p <^> spin at s p (p), Wp G TZ (8) 

simultaneously. If a site /x is orthogonal to a given root p, p ■ p = 0, then /x = s p (p). That 
is, affects only those spins at the sites p such that p ■ p ^ 0. In this way the action of the 
Hamiltonian (|j) in the model TZ is completely defined. 

The simplest example of the hierarchy is given by the anti- symmetric tensor represen- 
tation (AT) of A r corresponding to the Young diagram H. It consists of D = r(r + l)/2 
weights of the form 

g: AT = {e 3 +e k \ fc > j = 1, . . . , r + 1}, (9) 
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(1,2) (1,3) (1,4) (1,5) 

Figure 1: Structure of the second member of the A^p Heisenberg spin model hierarchy with 
nearest neighbour interactions. The spins on the sites which are connected by the same 
bonds are exchanged simultaneously. 

in which as before (|2]) we have omitted a constant vector 5 which is orthogonal to all the 
roots. For notational simplicity, we express the site ej + by a symmetric pair of integers 
(j, k) = (k,j) and the Pauli spin matrix on the site by (?(jk)- Obviously the sites {(j, fc)}'s 
form a two dimensional orthogonal lattice. A generic site (j, k) is connected to four nearest 
neighbouring sites: 

(j, k) (j - 1, k), (j + 1, k), (j, k-1), (j, k+1). (10) 

Those sites on the "boundary" 

boundary of £] = {(1, 2), . . . , (j, j + 1), . . . , (r, r + 1), (r + 1, 1)} (11) 

are connected to two other sites. The spin exchange operator V p for p = ej — exchanges 
the spins in the same row or column, i.e. (j,£) (k,£), W 7^ j, k, simultaneously: 

p = ej-e k , V p : e y + e t <$ e k + e t , i^j,k. (12) 

Thus for the anti-symmetric tensor (AT) model, the same Hamiltonian (|3]) no longer de- 
scribes the nearest neighbour interactions in the ordinary sense. As in the original Heisenberg 
model (^|) the spin exchange operator has a simple realisation in terms of the Pauli spin ma- 
trices: 

p = ej - e k =^ V p = [(1 + a m ■ a {k/) )/2]. (13) 

It should be emphasised that the two realisations fl5|) and (|13[) of the spin exchange operator 
V p (and any other realisations) satisfy the same commutation relations 

V P V*V P = V Sp (*), V 2 p = l, V^ P = V P , (14) 
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which are a direct consequence of the well-known commutation relations of the reflections 



SpSa^p ^s p {a)i $p 1) $ — p &p- 

Since the Hamiltonian commutes with the total spin operator 

S = \T,^ (16) 

z p,en 

the operators S 2 and S z define good quantum numbers. In the absence of the external 
magnetic field, the energy eigenstates consist of su(2) multiplets of spin S (2S + 1 fold 
degenerate) states as in the original Heisenberg spin chain. The vector (V), anti-symmetric 
tensor (AT) and the other higher rank anti-symmetric tensor representations of A r are special 
examples of the minimal representations which have played important roles in the Calogero- 
Moser Lax pairs [|13|, [14|, [15], [16], O] . The Calogero-Moser Lax pairs in these representations 



were first presented by D'Hoker and Phong ]nj. Minimal representations are characterised 
by the condition (see |T3, T5|) 



p v -fi = 0, ±1, WfieTZ, WpeA. (17) 

It is expected that the spin exchange models for these TZ's have some simpler features than 
others. 

In the rest of this paper we show the energy eigenvalues of the anti-symmetric tensor 
representation (AT) model for lower rank (r < 5) cases. We evaluate the Hamiltonian (|3|) 
in each S z = k, —r/2 < k < r/2 sector. In each sector, the Hamiltonian is a non- negative 
symmetric matrix with integer entries. In both V and AT models or any other models in 
the hierarchy, the highest eigenvalue is r + 1 which is the number of spin exchange operators 
V p in the Hamiltonian (|J). The lowest possible eigenvalue is — (r + 1). These two facts 
are based on the relation V"i — 1 ( |I4D implying that its eigenvalues are ±1. The highest 
eigenvalue r + 1 occurs in every S z sector but the lowest possible eigenvalue — (r + 1) appears 
only in smaller S z sector which have larger number of states. In particular, in the original 
Heisenberg (V) model, the lowest possible eigenvalue — (r + 1) never occurs, see (PH), (fZBD, 
(|27|) and ([3ll) . Once the lowest possible eigenvalue — (r + 1) is reached in an S z = k sector, 
then in all the \S Z \ < k sectors no new eigenvalues will appear. The energy eigenvalues of 
the vector representation (V) models are also listed for comparison. (The coupling constant 
J is put to unity for simplicity.) It is easy to see that the energy eigenvalues of the single 
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spin down (up) sector of the (V) model are given by 

2jn 



□ (I): ^ = r-l + 2cos 7 -^ TT , j = 0,l,...,r, 



1 



which are related to the eigenvalues of the Cartan matrix of and thus related to the mass 
spectrum of the afhne Toda field theory based on ||18||. All these eigenvalues appear in 



other sectors because of the su(2) invariance. It is also easy to see that for proper choices of 
the bases the Hamiltonian of the double spin down (up) sector of the (V) model is exactly 
the same as that of single spin down (up) sector of the (AT) model and so on: 

□ (II) =0(1), DUU) =Z U), (19) 

A 2 case In this case the vector rep. O (V) model and the anti-symmetric tensor rep. Q 
(AT) model are both 3-dimensional and these two models are equivalent from the spin model 
point of view. The total number of states is 8 = 2 3 and the energy eigenvalues together with 
(multiplicity) are 

3(4), 0(4). (20) 

The former consists of one S = 3/2 multiplet, whereas the zero energy states consist of two 
S — 1/2 states. 

A3 case The vector rep. □ is 4-dimensional whereas the anti-symmetric tensor rep. Q is 
6-dimensional. The spectrum of the □ (V) model is 

□ : 2 4 = 16 states 4(5), 2(7), 0(3), -2(1). (21) 

The maximal eigenvalue £ = 4 states are one S = 2 multiplet and the £ = 2 states consist 
of two S = 1 multiplets and one 5 = singlet. The £ = states are one 5 = 1 multiplet 
and the lowest energy £ = —2 state is a singlet. The spectrum of the Q (AT) model is 

[]: 2 6 = 64 states 4 (11), 2 (26), (12), -2 (14), -4 (1). (22) 

The maximal eigenvalue £ = 4 states are one S = 3, S = 1 multiplet each and a singlet. 
The £ = 2 states are three S = 2 and three S = 1 multiplets and two singlets. The zero 
energy £ = states are one S = 2 multiplet, two S = 1 multiplets and one singlet. The 
£ = —2 states have one S = 2 and three S = 1 multiplets. The lowest energy £ = —4 state 
is a singlet. 
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A A case The vector rep. □ is 5-dimensional whereas the anti-symmetric tensor rep. Q is 
10-dimensional. The spectrum of the Q (V) model is 

5± , 



32 states 5 (6), 



1 ± V5 (4) 



1 2 . 



(23) 



The maximal eigenvalue £ = 5 states are one S = 5/2 multiplet. The pair of energy levels 
£ = (5 ± v^5)/2 are both two S = 3/2 multiplets. Another pair £ = 1 ± V5 of states, one 
of which is the lowest energy, have each two S = 1/2 multiplets. The £ = 1 states are one 
S — 1/2 multiplet. The spectrum of the (AT) model is 

5± 



B = 



,10 



1024 states 5 (34 
-5± V$ 



;il2), liv^lie), -1±V5(76) 



32), ±V5 (36), 1 (58), (144), -1 (38) 



-5 (6). (24) 



This is a highly degenerate spectrum. The average multiplicity per energy level is 1024/15 ~ 
68 in contrast to 1024/139 ~ 7.3 in the vector rep. of Ay case (^TJ) which has the same number 
of sites as in this case. The highest £ = 5 states consist of 1 S = 5 (11-dimensional) , 1 
S = 3 (7-dimensional) ,25 = 2 (5-dimensional) and 25=1 (3-dimensional) multiplets, 
which we denote as follows: 



£ = 5 (34) : lxll, 1x7, 
Similarly we have the following degeneracy pattern: 



2x5, 2x3. 



(25) 



£ 


= (5±^5)/2 


(112) 


2x9, 


4x7, 


8x5, 


8x3, 


2x1, 


£ 


= i±VE 


(116) 


2x9, 


4x7, 


8x5, 


8x3, 


6x1, 


£ 


= -i±VE 


(76) 




4x7, 


4x5, 


8x3, 


4x1, 


£ 


= (-5± v / 5)/2 


(32) 






4x5, 


4x3, 




£ 


= ±y/E 


(36) 




1 x 7, 


3x5, 


4x3, 


2x1, 


£ 


= 1 


(58) 


1x9, 


2x7, 


4x5, 


4x3, 


3x1, 


£ 


= 


(144) 




4x7, 


12 x 5, 


16 x 3, 


8x1, 


£ 


= -1 


(38) 




2x7, 


2x5, 


4x3, 


2x1, 


£ 


= -5 


(6) 






1 x 5, 




lxl. 



(26) 



In the A$ models all the multiplicities are even numbers. 

A§ case The vector rep. □ is 6-dimensional whereas the anti-symmetric tensor rep. Q i s 
15-dimensional. The spectrum of the □ (V) model is 



64 states 



6(7 = 1x7), 5(10 = 2x5), 4(3=1x3), 3(10 = 2x5) 
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2(7 = 2x3+1x1), 1(6 = 2x3), 0(1), (27) 
(5±\/l7)/2 (6 = 2 x 3), l±V5(3 = lx3), l±Vl3(l). 

The spectrum of the Q (AT) model is 

Q : 2 15 = 32768 states 6 (156), 5 (776), 4 (1060), 3 (776), 2 (2756), 1 (4224) 

(1304), -1 (3600), -2 (1988), -3 (296), -4 (628), -5 (296), -6 (28), 
(5± Vl7)/2 (1200), (-5± vTr)/2 (576), 1 ± (600), -1±V5 (288), 
l±Vl3 (264), -l±Vl3 (240), ±c u ±c 2 , ±c 3 (1424), (28) 

in which 

Cj = (1 - lOcos[(0-2j'7r)/3])/3, j = 1,2,3, cos0 = 7/250. (29) 

These six eigenvalues always appear with the same multiplicities in each S z sector. This is 
a highly degenerate spectrum. The degeneracy pattern is as follows: 



£ 


total Mult. 16 


14 


12 


10 


8 


6 


4 


2 


6 


156 1 




1 


3 


4 


6 


6 


3 


5 


776 


2 


4 


10 


24 


34 


36 


30 


4 


1060 


1 


4 


13 


30 


50 


61 


42 


3 


776 


2 


4 


10 


24 


34 


36 


30 


2 


2756 


1 


8 


33 


76 


135 


168 


113 


1 


4224 


2 


12 


46 


118 


210 


256 


182 





1304 




3 


15 


34 


66 


85 


55 


-1 


3600 




6 


34 


100 


186 


234 


168 


-2 


1988 




2 


19 


52 


104 


136 


95 


-3 


296 






2 


8 


16 


20 


18 


-4 


628 






5 


17 


33 


45 


32 


-5 


296 






2 


8 


16 


20 


18 


-6 


28 








1 


2 


1 


2 


(5 ± y/V7)/2 


1200 


2 


6 


16 


34 


54 


64 


44 


(-5± Vl7)/2 


576 






1 


16 


30 


12 


30 


1±V5 


600 


1 


3 


8 


17 


27 


32 


22 


-1±V5 


288 






2 


8 


15 


21 


15 


i ± VT3 


264 




1 


4 


6 


13 


17 


9 


-1 ± 


240 




1 


3 


6 


12 


15 


9 


±Ci, ±Q2, ±C 3 


1424 




4 


11 


40 


74 


86 


64. 



The average multiplicity per energy level is 32768/31 ~ 1057. We believe this is a rule 
for the higher members of the spin model hierarchies. 
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Ag case The vector rep. □ is 10-dimensional which has the same number of sites as the 
anti-symmetric tensor rep. ^ °f ^4- There are 139 different energy levels among 1024 states. 
Here we list the spectrum (numerical values only) and (multiplicity) for comparison. 

□ : 2 10 = 1024 states 

10 (11), 9.618 (18), 9.345 (14), 9.173 (10), 9.118 (1), 9.099 (6), 9.064 (7), 8.982 (3), 
8.800 (10), 8.678 (14), 8.666 (5), 8.618 (18), 8.613 (10), 8.250 (6), 8.191 (14), 
8 (7), 7.978 (3), 7.963 (2), 7.934 (14), 7.855 (10), 7.800 (6), 7.594 (6), 7.421 (2), 
7.382 (18), 7.361 (10), 7.346 (10), 7.303 (10), 7.206 (10), 7.169 (6), 7.120 (14), 
6.893 (6), 6.730 (6), 6.695 (7), 6.677 (1), 6.653 (14), 6.540 (2), 6.501 (10), 
6.382 (18), 6.347 (6), 6.341 (3), 6.212 (1), 6.141 (6), 6.131 (6), 6.112 (10), 6.061 (10), 
6 (26), 5.902 (10), 5.899 (10), 5.708 (14), 5.627 (10), 5.618 (6), 5.459 (14), 5.364 (2), 
5.357 (14), 5.213 (10), 5.198 (6), 5.167 (10), 5.145 (10), 5.132 (6), 5.107 (6), 5.093 (2), 
5.018 (10), 4.929 (6), 4.902 (1), 4.848 (14), 4.666 (10), 4.648 (3), 4.646 (2), 4.618 (6), 
4.599 (5), 4.591 (6), 4.561 (10), 4.514 (10), 4.421 (14), 4.390 (10), 4.356 (6), 4.341 (2), 
4.232 (14), 4.171 (6), 4.163 (6), 4 (12), 3.939 (6), 3.697 (10), 3.566 (10), 3.563 (2), 
3.488 (2), 3.431 (2), 3.382 (6), 3.284 (10), 3.275 (6), 3.202 (6), 3.199 (10), 3.189 (6), 
3.091 (14), 3.027 (6), 3.025 (1), 2.964 (14), 2.869 (1), 2.747 (6), 2.694 (6), 2.687 (3), 
2.517 (10), 2.471 (10), 2.382 (6), 2.241 (7), 2.208 (10), 2.137 (10), 2.096 (3), 1.923 (2), 
1.916 (2), 1.904 (10), 1.696 (6), 1.687 (5), 1.646 (2), 1.565 (6), 1.523 (6), 1.454 (10), 
1.140 (6), 1.106 (6), 1.078 (6), 0.962 (1), 0.942 (6), 0.808 (1), 0.744 (10), 0.359 (10), 
0.332 (3), 0.122 (3), 0.025 (10), -0.029 (2), -0.535 (2), -0.561 (6), -0.773 (2), 
-0.782 (6), -0.952 (5), -1.492 (6), -2.087 (6), -2.541 (1), -3.184 (3), -4.031 (1). 

(31) 

Some remarks and comments are in order. 

One might feel that the simultaneous exchange of the spins is artificial. Rather it is quite 
natural within the context of the conventional spin exchange models, if we consider any 
member of the conserved quantities as a new Hamiltonian but keeping the sites of the spin 
the same as in the original model. As shown in the previous paper |l|] the conserved quantities 
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consist of the linear combination of ordered products of the spin exchange operators. Each 
product term of the conserved quantity specifies simultaneous exchange of spins. The only 
exception is the original spin exchange models based on A r root systems and the vector 
representation V. For the other root systems, even for the classical ones B r , C r or D r , the 
single spin exchange operator implies a multiple exchange as shown explicitly in (3.21) of 
I- 

It would be interesting to determine the energy eigenvalues of the Hamiltonian ([3|) for a 
different choice of the set of affine simple roots, for example D^: 

n (1) = {ex - e 2 , e 2 - e 3 , . . . , e r ^ - e r , e r _i + e r , -(e x + e 2 )}, (32) 

and for the D r vector representations or (anti-) spinor representations. The former would 
give a .D-series analog of the nearest neighbour Heisenberg model. 

Another well-known method for spin models is the use of Yang-Baxter equations [JTJJ. It 
would be interesting to compare the fusion procedure for the i?-matrices with the present 
method of generating hierarchies. 
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